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We consider sample of inclusive 7 + 3 jet events collected by the DO experiment. The double parton fraction 
(/dp) an <i effective cross section cr e ff , a process-independent scale parameter related to the parton density inside 
the nucleon, are measured in three intervals of the second (ordered in pt) jet transverse momentum p J ^ t2 within 
the 15 < J^y* 2 < 30 GeV range. Also we measured cross sections as a function of the angle in the plane 
transverse to the beam direction between the transverse momentum (px) of the 7+leading jet system and px of 
the other jet for 7 + 2 jet, or pt sum of the two other jets for 7 + 3 jet events. The results are compared to 
different models of multiple parton interactions (MPI) in the PYTHIA and SHERPA Monte Carlo (MC) generators. 



1. Introduction 

Many features of high energy inelastic hadron collisions depend directly on the parton structure of hadrons. 
The inelastic scattering of nucleons occurs mainly through a single parton-parton interaction but the contribu- 
tion from double (or multiple) parton collisions can be significant. Information about DP rates is needed for 
understanding of nature of MPI events and correct estimating background to many rare processes, especially 
with multi-jet final state. 



2. Double parton interactions in 7+3 jet events 

The cross section of DP events production directly proportional to cross sections of two processes A and B 
and should be normalized by some scaling parameter a e s in cross section's units. 

o DP = . (1) 

Ceff 

In general sense, c e g is a factor which characterizes a size of the effective interaction region. 

We have used a sample of 7 + 3 jets events collected by the DO experiment with an integrated luminosity of 
about 1 fb — 1 . The DO detector is a general purpose detector described in The events should pass triggers 
based on the identification of high px cluster in the EM calorimeter with loose shower shape requirements for 
photons. Jets are reconstructed using the DO Run II iterative midpoint cone algorithm 0| with a cone size 0.7. 
Each event must contain at least one 7 in the rapidity region \y\ < 1.0 or 1.5 < \y\ < 2.5 and at least three 
jets with \y\ < 3.0. Events are selected with 7 transverse momentum 60 < p T < 80 GeV, leading (in pr) jet 
Pt > 25 GeV, while the next-to-leading (second) and third jets must have px > 15 GeV. The DP fractions and 
cr e ff are determined in three p^ t2 bins: 15-20, 20-25, and 25-30 GeV. 

We use rates of double interactions (DI) from two separate pp collisions and DP from a single pp collision to 
extract <7 e ff from their ratio 3]. The DI events differ from the DP events by the fact that the second parton 
scattering happens at a separate pp collision vertex. Data events with a single pp collision vertex, which compose 
the sample of DP candidates, are selected separately from events with two vertices which compose a sample of 
DI candidates. A distinctive feature of the DP events is a presence of two independent parton-parton scatterings 
within the same pp collision. We define the variable sensitive to the kinematics of the DP events: 

AS = A0 (f T ' jet \ 4 et2je * 3 ) , (2) 

where A(f> is an azimuthal angle between the px vectors of the total transverse momenta of the two two-body 
systems, p T ^ etx and p^f t2je * 3 , in 7 + 3 jet events. This angle is schematically shown in Fig. EJ 

We consider the data-driven method to extract the DP fractions /dp- Specifically, we consider data in two 
adjacent px intervals of the second jet. The distribution for AS variable in data can be expressed as a sum 
of signal and background distributions. If we known properties of data and DP model, the only unknown 
parameter is the fraction of DP events in one p^ t2 bin. It is obtained from a minimization. The found /pjp 
values with total uncertainties are 0.466 ±0.041 for 15 < p^ t2 < 20 GeV, 0.334 + 0.023 for 20 < p^ 2 < 25 GeV, 
and 0.235 ± 0.027 for 25 < p^, t2 < 30 GeV. They are shown in Fig. [5] (three sets of the points correspond to 
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Figure 1: A possible orientation of photon and jets transverse momenta vectors in 7 + 3 jet events. Vectors and 
are the pr imbalance vectors of 7+jet and jet-jet pairs. The figure illustrates a general case for the production of 7+3 
jets +X events. 



three possible definitions for the A5* variable [3] ) . The values of a e s are shown in Fig. [3] The main systematic 
uncertainty are caused by determinations of the DI and DP fractions giving a total systematic uncertainty of 
(20.5 — 32.2)%. The obtained <7 e ff values in different p!^ t2 bins agree with each other within their uncertainties 
and highly uncorrelated, and are used to calculate the average value: 

cr^ c = 16.4 ± 0.3(stat) ± 2.3(syst) mb. (3) 

This average value is in the range of those found in previous measurements 0-0] performed at different energy 
scales of parton interactions. 
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Figure 2: Fractions of 7 + 3 jetevents with double parton interactions in the three p^ t2 intervals. 
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Figure 3: Effective cross section <r G fi (mb) measured in the three J3y t2 intervals. 



3. Azimuthal decorrelations and multiple parton interactions 

in 7 + 2 jet and 7 + 3 jet events in pp collisions 

Samples of 7 + 2(3) jets events with the same cuts as [3| are considered. The next modifications are applied: 
each event must contain at least one 7 in the pseudorapidity region \y\ < 1.0 or 1.5 < \y\ < 2.5 and at least 
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two (or three) jets with |y| < 3.5. Events are selected with 7 transverse momentum 50 < jP T < 90 GeV, leading 
jet pt > 30 GeV, and the second jet pr > 15 GeV. If there is a third jet with pt > 15 GeV that passes the 
selection criteria, the event is also considered for the 7 + 3 jet analysis. By measuring differential cross sections 
vs. the azimuthal angles in 7 + 2(3) jet events we can better tune MPI models in events with high px jets. 
We present the four measurements of normalized differential cross sections, AS in a sing le bin 15 < p!^ < 30 
for 7 + 3 jet events (see Fig. [IJ, and Atfi in three pl^ t2 bins, 15-20, 20-25, and 25-30 GeV, for 7 + 2 jet events. 
The A(j) is an angle between the pt vector obtained by pairing the 7 and the leading jet pr vectors and the 
second jet px vector Q). It is shown in Fig. 2) We consider a few MPI models and two models without MPI 
simulated by PYTHIA and SHERPA MC generators. Figure [5] shows the measured cross section for the two angular 
variables AS" (left plot) and A(j> (right plot). The data have a good sensitivity to the various MPI models, which 
predictions vary significantly and differ from each other by up to a factor 2 at small AS and Ac/), i.e. in the 
region where the relative DP contributi 




Figure 4: A possible orientation of photon and jets transverse momenta vectors in 7 + 2 jet events. 
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Figure 5: Left: Normalized differential cross section in the 7 + 3-jet events, (l/o" 7 3j)a" 7 3j /dAS , in data compared to 
MC models and the ratio of data over theory, only for models including MPI, in the range 15 < p^ t2 < 30 GeV. Right: 
Normalized differential cross section in 7 + 2-jet events, (l/rj 7 2j)o"72j/dA0, in data compared to MC models and the 
ratio of data over theory, only for models including MPI, in the range 15 < p J T et2 < 20 GeV. 



From these two plots we may conclude that: (a) a large difference between single parton-parton interaction 
(SP) models and data confirms a presence of DP events in the data sample; (b) the data favor the predictions 
of the MPI models with P0, SO and Sherpa MPI tunes with py-ordered showers; (c) the predictions from tune 
A and DW MPI models are disfavored. It is important that our preferable choice of MPI models is stable for 
all our measurements. 

In 7 + 2 jet events in which the second jet is produced in the additional independent parton interaction, the 
Acf> distribution should be flat. Using this fact and also SP prediction for A<fi we can get the DP fractions from 
a fit to data. The distributions in data, SP, and DP models, as well as a sum of the SP and DP distributions, 
weighted with their respective fractions for 15 < p^ t2 < 20 GeV, are shown in the left plot of Figure [6] The DP 
fractions in the 7 + 2 jet samples decrease in the bins of p^ t2 as (11.6 ± 1.0)% for 15 — 20 GeV, (5.0 ± 1.2)% for 
20 — 25 GeV, and (2.2 ± 0.8)% for 25 — 30 GeV. To determine the fractions as a function of Acf>, we perform a 
fit in the different A<f> regions by excluding the bins at high A</>. We find that they grow significantly towards 
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the smaller angles and are higher for smaller p^, t2 (right plot of Figure . 
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Figure 6: Left: the A(f> distribution in data, SP, and DP models, and the sum of the SP and DP contributions weighted 
with their fractions for 15 < p J ^ t2 < 20 GeV. Right: the fractions of DP events with total uncertainties in 7 + 2 jetfinal 
state as a function of the upper limit on A(j> for the three p J ^ t2 intervals. 

We also estimate the fraction of 7 + 3 jet events from triple parton interactions (TP) in data as a function of 
p^ t2 ■ In 7 + 3 jet TP events, the three jets come from three different parton interactions, one 7+ jet and two 
dijet final states. In each of the two dijet events, one of the jets is either not reconstructed or below the 15 GeV 
Pt selection threshold. The fractions of TP events in the 7 + 3 jet data have been estimated and are shown in 
Fig. As we see, they vary in the p^ t2 bins as (5.5 ± 1.1)% for 15 - 20 GeV, (2.1 ± 0.6)% for 20 - 25 GeV, 
and (0.9 ± 0.3)% for 25 - 30 GeV. 
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Figure 7: Fractions of 7 + 3 jetevents with triple parton interactions in the three p J y t2 intervals. 
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